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ABSTRACT 

We present clear spectroscopic detections of molecular hydrogen in the outer filaments of the 
Ha nebula surrounding the central galaxy of the Perseus cluster, NGC 1275. This implies the 
presence of warm molecular gas clouds at projected distances up to 24 kpc from the nucleus of 
the host galaxy, embedded in the hot intracluster medium. The emission-line intensity ratios 
reveal that the H2 emission is predominately thermally excited with excitation temperatures 
of 1600— 2200 K, suggesting a lack of pressure balance between the molecular component 
and its surrounding medium. Excitation by stellar UV or by the central AGN is shown to be 
unlikely, whilst thermal excitation by X-rays or conduction from the ICM or shocks are all 
possibilities. Evidence for a non-thermal component is found in the spectra of some regions 
based on a low ortho-to-para ratio. 
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1 INTRODUCTION 

Luminous optical and UV line-emitting filamentary nebulae like 
that of NGC 1275 are commonly found surrounding central cluster 
galaxies where the X-ray emitting intracluster medium (ICM) has 
a short radiative cooling time ICrawford et all 19991 and references 
therein). Despite many studies, the origin and excitation mecha- 
nism of these filaments are uncertain, with cooling of the intra- 
cluster gas, star formation, conduction and shock heating all hav- 
ing been suggest ed and contested (e.g. Ijohnstone & Fabianlll988[ 
ISabra et all200ch . 

The large Ha nebulousity surrounding NGC 1275 has been 
long known 1 Minkow ski! 1 19571 lLvndsMl970j) and has been stud- 
ied in a variety of wavelengths. Dee p Hq and X-ray images ob- 
tained by the WIYN telescope JConselice etai]|200lh and Chan- 
dra ( Fabia n et al]|2003l) show spectacular, extended Ha filamen- 
tary structures stretching over 50 kpc from the central galaxy, with 
corresponding soft X-ray counterparts, deeply embedded in the hot 
ICM. 

NGC 1275 lies at the centre of a cluster with strong centrally 
peaked X-ray emission where th e X-ray luminosity implies a mass 
deposition rate of 300M e yr _1 I Allen et af]|200lh . However, after 
extensive searching, the massive amounts of gas expected to accu- 
mulate cannot be accounted for as stars or line-emitting gas near 
the nucleus. 

Chandra spectra show little evidence of gas cooling below 
a temperat ure about one third of that in the bulk of the cluster 
( Schm idt et ai]|2002t lianders et alJl2004l) . Therefore most of the 
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gas either cools further in a non-radiative manner, requiring en- 
ergy to be shared between the hot X-ray emitting and cooler gas 
by mixing or conduction, or is heated by a distributed heat source. 
A number of theories including heating by the central active galac- 
tic nu cleus (AGN) of NGC 1275 have been suggested (see Fabian 
2004 and references therein). 



A reservoir of ~4x 10 6 M faeckman et"ai]l 19891) of 10 4 - 
10 5 K gas lies in the optical and UV line-emittin g nebula. Molecu - 
lar gas appears common in such regions jjaffe & Breme3ll997h . 
Recent observations show that there is a compone nt of warm 
(1000-4000 K) molecular hydrogen gas iDonahue et ail l2000t 
Edgeetal. 2002) with up t o ~1 x 10 10 Ma of cool mo lecular gas 
lEdgell200ll) iKrabbe et ail feOOOl) and lDonahue et alj <2000h de- 
duced that the molecular hydrogen component in the NGC 1275 
nucle ar region is lik ely to be excited by the central AGN, whilst 
IWilman et alj 120021) argue that the H2 emission is thermal and 
mostly likely to be heated by UV radiation. The outer filaments 
have a fairly constant Ha surface brightness extending over many 
tens of kpc, therefore a central excitation source such as the AGN is 
unlikely to be of great importance. Whilst massive stellar clusters 
are seen in Ha, they do not spatially correspond with all the fila- 
ments and bright knots. Thus the excitation of the outer filaments 
may be different to that of the nuclear region. 



The redshift of NGC 1275 is 0.0176, 
Hq=70 kms~ 1 Mpc~ 1 , gives 1 kpc~2.7 arcsec. 



which using 
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Target Slit centre co-ordinates Position Nodding Total integration 

RA Dec (J2000) angle (arcsec) time on target (minutes) 

Slitl(NW) 03 1946.78 +41 31 45.28 -125° 0,-10 36 

Slit2(SW) 03 19 41.95 +41 29 58.61 -76° -30,-20 48 

Slit 3(E) 03 19 49.31 +41 31 50.52 -4° -35,-80 54 



Table 1. Log of Observations. Nodding 1 direction and position angle as explained in CGS4 literature. 



2 OBSERVATIONS 

The spectra presented in this paper were taken with the CGS4 spec- 
trograph on the United Kingdom Infrared Telescope on 2004 Jan- 
uary 17, 18 and 19 for a total integration time of 15 hours. The 
40 lines mm~ 1 grating was used with the 300 mm focal length cam- 
era, giving a spatial scale of 0.61 arcsec per pixel. The B2 filter cen- 
tred at 2.02 |j.m was used with a two pixel wide slit giving a spectral 
resolution of 600—790 kms~ 1 . The observations were taken in ND- 
STARE mode, using a object-sky-sky-object nodding pattern. The 
seeing was better than one arcsec on all three nights and there was 
little to no cloud cover. 

The telescope was directed to a bright reference source near 
the target position. The slit was then shifted to a targeted bright 
object (co-ordinates given in Table0 and orientated to the correct 
position angle. For confirmation of the exact spatial position of the 
slit we included one or two bright objects such as a star or galaxy 
along the length of the slit. TableQgives the co-ordinates, position 
angle and integration time for all observations. Guided by the Ha 
map of Conseli ce et alJ 1200 ll) . the slit positions were chosen to 
include some of the brightest knots and filaments furthest from the 
nucleus. 

Most sky emission features were removed automatically by 
the nodding pattern. An atmospheric, temporally variable CO2 line 
near 2.01 fim and a feature at 2.15 \im however remain in the 
spectrum from the Slit 1 and 2 position. The target observations 
were preceded and followed by a series of calibration observa- 
tions including flat-fields, dark frames, arcs and a set of standard 
F and G type stars. These standards were nearly featureless stars 
of known magnitude. In each of the standard stars a Bry stellar 
feature at 2.166 fim was removed by linear interpolation. The data 
was reduced with ORAC-DR Version V4.0-1 available through the 
UKIRT software website. Data manipulation and extraction was 
done with Starlink packages, QDP iTennant 1990) and IDL Ver- 
sion 6.0. 



3 RESULTS 

Clear detections of the strongest molecular hydrogen lines: 
n=l—0 S(l) and ?;=1— S(3) were found in seven separate regions 
of the outer filaments, marked and labelled in Fig. Re gion A 
covers an area of very strong Ha emission. Shields & FiliDDenko 
1 1990) discovered and studied a massive, young stellar cluster in 
this region. It lies 19.4 arcsec (a projected distance of 7.2 kpc) to the 
East of the central galaxy, marked as a bright knot in Ha within the 
Northern bright filament that is part of array of radial filaments that 
point East from the galaxy. Regions B, C and D come from this ar- 
ray of radial filaments at distances 18.0 arcsec (6.7 kpc), 17.7 arcsec 
(6.6 kpc), 17.9 arcsec (6.6 kpc) from the central galaxy respectively. 
There is no clear distinction separating regions C and D in the slit 
image, although spectral analysis revealed they were most likely 
separate populations of molecular hydrogen. The 'horseshoe knot' 



Region 

RA Dec (J2000) 

size of extracted region 



Line 



Surface Brightness 



Horseshoe knot 


o=l 


-0S(3) 


20.7 +jjf 


03 19 45.22 +41 31 32.61 


v=l 


-0 S(2) 


6 6 +2 S 
DD -2.5 


2.23 arcsec 2 


o=l 


-0S(1) 


18.1 




v=l 


-0 S(0) 


5 2 +3[6 


SW1 knot 


Pa 




11 5+ 7 ' 2 


03 19 45.69 +41 29 46.99 


v=l 


-0 S(4) 


7 5 + 2 - 9 
-2.7 


4.47 arcsec 2 


o=l 


-0 S(3) 


13.1 




o=l 


-0 S(2) 


2.6 ±1.4 




v=l 


-0S(1) 


15.31" 




v=l 


-0 S(0) 


4-9 t\t 




o=l 


-0Q(1) 


14-7 111 



SW2 knot 
03 19 44.97 
2.98 arcsec 2 



+41 29 48.92 



o=l-0S(3) 
ti=l-0S(l) 



Region A 

03 19 49.71 +41 30 51.11 
2.23 arcsec 2 



Pee (Narrow) 
Pa(Broad) 
v=l-0 S(3) 
o=l-0 S(l) 
B17 



Region D 
03 19 49.77 
2.23 arcsec 2 



+41 30 41.28 



d=1-0S(3) 
«=1-0S(1) 



4.6; 



,+2.7 
1 -2.2 

5 6 + 21 
3.0 _ L9 



119.1 ±35.6 
38.7 ±35.6 

46+13 
t-o _i. 4 

6 5 + 17 

-2.6 



Region B 


Pa 


4i ,+10.9 

^' -9.0 


03 19 49.75 +41 30 45.19 


o=l-0 S(3) 


23.5 tH 


2.23 arcsec 2 


o=l-0 S(2) 


12 5 +2 - 6 




o=l-0 S(l) 
o=l-0 S(0) 


24.5 + " 
9-8 IS 




o=l-0 Q(l) 


17.21" 




o=l-0 Q(3) 


18.6 



Region C 


Pa 




25-5 l«;t 


03 19 49.77 +41 30 43.24 


o=l 


-0S(3) 


16.8 ±1.9 


2.98 arcsec 2 


o=l 


-0 S(2) 


7 3 +L9 




o=l 


-0S(1) 


10.5 ±1.6 




o=l 


-0 S(0) 


4 6+17 
-1.6 




o=l 


-0Q(1) 


10 2 +4 ' 4 



10.3 



+3.0 

' -2.5 

5 1 +u 
3 - 1 -0.4 



Table 2. Line surface brightness in 10 erg s cm arcsec for all re- 
gions (as marked in Fig.Q in which the strongest ro- vibrational H9 lines, 
o=l—0 S(l) and o=l— S(3) are clearly detected. Errors are at 1 o~ level. 



lies at the tip of a horseshoe- shaped feature at the end of a 74 arcsec 
lo ng NW radial filament t hat appears to turn back on itself, noted 
by IConselice et alJ 1200 ll) (box 1 1 in their Fig. 5). It has a pro- 
ject ed distance of 22 kpc away from the galaxy center. SW1 (box 
2 in lConselice et all200ll Fig. 5) and SW2 are bright knots that lie 
at the end of radial filaments (61.6 arcsec ~23 kpc and 64.6 arcsec 
~ 24 kpc from the galaxy respectively) pointing SSW. Our K-band 
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Figure 1. Left: WIYN telescope Ha image of the line-emitting filamentary structure surrounding NGC1275 from lConselice et all 1200 lh . Regions where H2 
ro- vibrational lines are detected have been marked by circles and labelled. Right: Detail of the Eastern slit region. 
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Figure 2. UKIRT K-band smoothed spectra of the outer filament regions marked in Fig.Qwith the ro-vibrational H2 and atomic hydrogen lines labelled. Line 
labels have redshift z=0.0172. Errors are from Poisson Statistics. Note that the scale in Region A does not show full extent of Pa. 
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spectra for all regions are presented in Fig. [2] with line surface 
brightness given in Table|2| Since the structures observed may be 
smaller than the slit size, these intensities may be underestimated, 
hence they could be considered lower limits. The spectra have not 
had any stellar continuum removed. The SW regions and the horse- 
shoe knot show no evidence for a continuum, whereas regions A, B, 
C and D all have a continuum due to the galactic stellar emission. 

In principle the wavelength range should allow an esti- 
mation of the extinction using the Pa and Bry lines or from 
the n=l-0 Q(3) and v=l— S(l) lines. However the Pa and 
1=1—0 Q(3) emission lines pass through complex parts of the at- 
mospheric window and the Bry emission is too faint to have a sig- 
nificant detection in all but one region. Therefore no reddening has 
been taken into account in this work, although we note any conse- 
quence it may have in our discussion of the data. 




19 191 1.92 1.9 1.91 1.92 

Wavelength (^m) Wavelength (^im) 



Figure 3. Pa line emission from region A, the young stellar cluster (left) 
and region B (right), both showing non-Gaussian profiles. 



4 ANALYSIS AND DISCUSSION 
4.1 Velocity Structure 

The Pa emission from regions A and B both exhibit non-Gaussian 
profiles, shown in Fig. [5] This profile could not occur through line 
blending as the Pa line is distinct from the other molecular and 
atomic hydrogen lines. These profiles may be the result of either 
extra redshifted material or absorption of the blueshifted emission, 
either intrinsically or by the Earth's atmosphere (this line is located 
in a complex part of the atmospheric window). Neither the molec- 
ular hydrogen lines nor the other atomic lines from our spectra ex- 
hibit these profiles, although given the high Pa/H2 line ratio of 
region A, it is unlikely that a low flux red wing would be visible in 
the molecular emission above the noise level. 

The Ha emi ssion from region A has been studied by 
IShieldsetalJ<1990h who find the Ha profile shows a well defined 
narrow peak at z=0.0168 and a large bump of emission at a larger 
redshift. They interpret the lower redshift, narrow emission as orig- 
inating from the stellar cluster, and the broad emission originat- 
ing from overlapping filamentary material. From the [Ol]/Ha ratio 
they estimated approximately 30 per cent of the total Ha luminos- 
ity in this region originates from the filament. The Pa emission 
profile from region A can be better fitted by two Gaussians, a nar- 
row Gaussian at z~0.0165 containing 75 per cent of the total line 
emission and a broad Gaussian containing 25 per cent of the line 
emission with z~0.018, rather than a single Gaussian. Although 
ou r slit may not have covered the exact same region investigated 
by Shi elds et alJ ll990h . the similarities in line profile suggest the 
same interpretation. 

Region B can also be fitted by two Gaussians; a broad domi- 
nating component at z~0.0169 and a narrow component at a much 
larger redshift. The large errors and lack of other corroborating evi- 
dence makes it difficult to apply an intrinsic emission interpretation 
to this region as opposed to absorption by the atmosphere. There- 
fore a single Gaussian profile is fitted to the Pa emission. 

None of the lines were resolved spectrally but Fig. [4] shows 
that there is a velocity shift relative to the nucleus of the strong 
H 2 lines: «=l-0 S(l), «=l-0 S(2), w=l-0 S(3) and Pa emission. 
This figure shows the molecular and atomic emission from regions 
B, C and the SW1 knot appear to have similar velocity shifts, there- 
fore it is likely they share an origin. The Bry emission from region 
A has a similar velocity shift to the narrow Pa emission, therefore 
most likely to have originated from the same region, the stellar clus- 
ter. The molecular emission from region A does not share the same 
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Figure 4. Velocity shift of molecular and atomic hydrogen lines relative to 
the nucleus, measured from a combined fit of the v=l— S(l), v=l— S(2) 
and v=l—0 S(3) lines (crosses), Pa lines (diamonds), Bry (square). HS is 
the horseshoe knot. The velocity shift of both the narrow Pa emission la- 
belled A(SC) and the broad emission labelled A(fil) are shown. Velocity 
resolution is approximately 700 kms , errors reflect 1 o" level on line posi - 
tion. Nuclear molecular hydrogen line positions from lKrabbe et al . I.2OO0); 
nuclear Pa (1 .8756 \xm) and Br/ (2. 1661 yeca) line positions chosen to con- 
form with redshift of z=0.0172 as measured bv lKrabbeet"ai]|200dl 



velocity shift as either the stellar cluster or the filament. It is possi- 
ble the molecular hydrogen emission is a superposition of emission 
from both the filament and the stellar cluster, or the molecular emis- 
sion may not share the same origin as the atomic emission. If the 
Pa emission was fit with a single profile, the velocity shifts of the 
Pa and the H2 would still not agree. This figure suggests that the 
outer filaments observed are systematically blueshifted relative to 
the nucleus. 



4.2 H2 Excitation 

H2 excitation may occur through collisions and therefore ex- 
hibit thermal emission, or from radiative transfers producing non- 
thermal emission. These mechanisms can be distinguished since 
they preferentially populate different levels producing different 
spectra, but as they often occur together the observed spec- 
tra can exhibit both thermal and non-thermal emission. Colli- 
sional de-excitation dominates when the gas has a total den- 
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Figure 5. Diagram constraining likely excitation mechanisms. The four 
3o~ upper limits indicate observed line ratios from the brightest fila- 
ments, errors on the v=l— S(1)/d=1— S(3) ratio are lo~ level. Boxes 
highlight reg ions occupied by the non-thermal UV excitation models of 
iBlack & van DishoectJ fl987) and models o f thermal processes near the 
dotted thermal curve including UV excitation {Sternberg & Dalgarnol 19891) 
and X-ray heating models i Lepp & McCrav 1983 and Draine & Woods 
Il99fil . The x-axis gives an indication of the level of thermal to non-thermal 
excitation occurring. The upper limits are consistent with a thermal excita- 
tion mechanism, but there also may be a non-thermal component. 



sity n T = n H +n H , > 10 4 cm" 3 JSternberg & Dalgarnol Il989t 
iMandv & Martmlll993h for T~10 3 K. The gas may be heated by 
shocks, conduction from the ICM, X-rays or UV photons and the 
H2 energy levels populated through collisional transitions. When 
this occurs the ortho-to-para ratio is three and the gas is in local 
thermodynamic equilibrium (LTE). Below the critical density, the 
probability of the H2 undergoing a transition from an upper to a 
lower energy level via a radiative process is greater than for colli- 
sional de-excitation. The H2 molecules can be excited to the upper 
energy levels through non-thermal excitation mechanisms such as 
X-rays and UV photons in the Lyman-Werner band (912— 1108A) 
or through collisions with suprathermal electrons caused by X-ray 
ionization. 

Non-thermal excitation mechanisms readily excites the v = 
2 and higher vibrational states, whereby collisional transitions 
preferentially de-excite the v = 2 level in favour of the v = 1 
level, so the v=2—l S(l)/ v=l— S(l) ratio is an indicator of 
the relative level of thermal to non-thermal processes. Reddening 
does not have a large effect in this part of the spectrum, there- 
fore it can be ignored without much consequence when consider- 
ing H2 line ratios. Models predict a v=2—l S(l)/«=1— S(l) ra- 
tio greater than 0. 53 for pure fluorescence at rij ^10 5 cm~ 3 and 
low temperatures I Draine & Bertoldi 1996), whilst shock models 
of Hollenb ach & McKeel ll989) exhibit ratios between 0.04—0.46. 
3o~ upper limits of the v=2— 1 S(l) emission from the four regions 
exhibiting the strongest H2 lines and best signal-to-noise ratio (re- 
gion A, region B, SW2, horseshoe knot) are plotted in an excitation 
diagram (Fig.[3J. All lines used in these ratios are ortho, so this di- 
agram is independent of an uncertain ortho-to-para ratio. Fig.[5]is 
consistent with predominately thermal (collisional) excitation oc- 
curring but does not rule out a non-thermal component. 

If collisional de-excitation is assumed to dominate, the H2 will 
be in local thermal equilibrium (LTE). The energy levels of H2 will 
be populated in a Boltzmann distribution and from this an excita- 



N(v,J) 
gJ 



-E{v,J)/k b T a 



where E(v,J) is the upper energy of the (v,J) transition, a is a 
constant, kj, is Boltzmann's constant and gy is the statistical weight. 
N(v,J) is the transition column density, which, for optically-thin 
emission is given by 

N(v,J) = ^L 
' A ul / ? v 

where / is the observed intensity of the line, v is the rest fre- 
quency and A u i is the Einstein co-efficient taken from lTurner et alJ 
1 1977). In all calculations an ortho-to-para ratio of 3:1 has been 
assumed. If the H2 is only thermally excited, the excitation temper- 
ature is the kinetic temperature. However, if non-thermal excita- 
tion is also taking place, the excitation temperature is not the actual 
temperature of the gas, but characterises the level populations that 
have arisen through competing thermal and non-thermal processes 
IPuxlev et alfcOOfj) . 

The LTE plots (Fig.|6j show the horseshoe knot spectrum fits 
a model of thermalised gas with an ortho-to-para ratio of three. The 
LTE plots of regions B and C show a clear misalignment of the 
ortho (odd J) and para (even J) lines which indicate a lower ortho- 
to-para ratio. The v=\— S(4) line (and possibly v=l— S(2)) are 
much narrower than the other molecular lines from the SW1 region, 
possibly affected by atmospheric sky subtraction (see Fig. 2). The 
S(0) line appears to agree with a model of thermalised gas with an 
ortho-to-para ratio of three, whilst the S(2) and S(4) lines do not. 

Molecular hydrogen is formed with an ortho-to-para ratio of 
three. No radiative process can induce conversions between the or- 
tho and para lines, hence the ratio is conserved. There are three 
known methods in which ortho-to-para conversions can occur: by 
collisions with protons in cold (T<100K) gas, or collisions with 
hydrogen atoms in warmer gas (T<300K). Finally, if the hydro- 
gen molecule stays on a grain long enough it can evaporate with 
an ortho-to-para ratio appropriate to the grain temperature. In re- 
gions where shock excitation is the primary excitation mechanism, 
the rapid collisio ns cause the orth o-to-para ratio to take on a value 
close to three ISmifh et alJ ll997). However photodissociation re- 
gions (PD Rs) often e xhibit lower ortho-to-para ratios, typically 
1.5-2.2. ISternberg & Neufekl < 19991) provides models of far-UV 
(FUV) pumping in PDRs with optically thick H2 (~10 14 cm~ 2 for 
UV) which explain how the observed ortho-to-para ratio will tend 
to 1.7 whilst the true ortho-to-para ratio remains at three due to 
lower rates of FUV pumping of molecular hydrogen in the optically 
thicker ortho states compared to the para states. In these models, the 
gas density must remain below 5x 10 4 cm~ 3 and be relatively cool 
(T< 1000K) in order for collisional excitation to be negligible rela- 
tive to FUV pumping, and collisional de-excitation to be negligible 
compared to spontaneous radiative decay, otherwise the ortho-to- 
para ratio would tend t o the LTE value of three. 

We follow Draine & Bertoldi 1 1996) definition of the ortho- 
to-para ratio, y 



This ratio can be determine d fro m a single pair of one ortho and 
one para line JPuxlev et all200fj) through 

Io A p (2J p + 1)A cAEp j kbTa 
' l p A (2J + \)X p 

where X is the wavelength, J is the rotational quantum number 
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Figure 6. Local thermodynamic equilibrium plots, where the inverse slope of the best fit line gives T ex . Errors on data points are at 1 o~ level, and error on T ex 
is Iff error from weighted best fit. The labels refer to the d=1— S and Q transitions. v=l— S(l) and v=l—0 Q(3) transitions have been offset in region B for 
clarity. 



of the upper state of the transition and AE po is the energy differ- 
ence between the upper energy level of the para and ortho line. 
Ortho/para lines are denoted by a o/p subscript. 

T ex and ortho- to para ratios are given in Table|3|for regions B, 
C, SW1, and the horseshoe knot, where three or more Irb lines are 
detected. The excitation temperature of the molecular h ydrogen gas 
in the centre of NGC 1275 is given bv lWilman et all 120021) to be 
1620 K, whereas Krab be et alJfeOOOh measure it to be 1450±250 K 
(from the v = 1 transitions). T ex of the filaments are consistent with 
these values for the gas in the core of the nebula, although tend to 
be slightly higher. 

Although there are large errors, the ortho-to-para ratios sug- 
gests that regions B and C may also be non-thermally excited, pos- 
sibly by FUV pumping in a PDR. 

Assuming the upper levels are also populated in a Boltzmann 
distribution at the same excitation temperature, the total H2 column 
density is given by 



^Total 



N{v,J)Z(T) 

-E(vJ)/k b T a 



gje 



where Z(T) is the partition function. This expression assumes all 
the molecular hydrogen is collisionally excited at the same tem- 
perature and all the emission is visible. There are many scenarios 
in which this may no t be the case, such as the X-ray dissociated 
regions modelled bv lMalonev et alJJl996» . In addition the line sur- 
face brightness listed in Tableware only lower limits, therefore 
these column densities can be considered lower limits also. 

Column densities are listed for each region with three or more 
strong H2 lines in Table|5| If T ex of the upper levels is greater than 



in the lower levels, the column density will decrease from these 
values. All column densities are low (2—14 x 10 15 cm~ 2 ) which 
has implications for the structure of the filaments. For the gas to 
be thermalised the total density must be greater than 10 4 cm~ 3 . 
The Ha filaments are estimated to have a density of ~10 2 cm -3 
assuming they are in pressure equilibrium with the surrounding 
X-ray emitting gas. The filaments must thus have an undetected 
denser and cooler (1600-2200 K) atomic component that mixes and 
shares energy with the molecular hydrogen. Alternatively, the H2 
exists in clumps, much smaller than a parsec in diameter, strung 
through the Ha filaments. This has parallels with suggestions of 
Fabia n~et alJ 120031) that the atomic hydrogen has a similar struc- 
ture, on the basis of the Ha emission and pressure balance. Con- 
sistently Wilman et al. 1 2002) argue that the H2 emission near the 
nucleus of NGC 1275 must come from a population of dense self- 
gravitating clouds based on H2 line ratios. 



4.3 Comparison of Molecular and Atomic Hydrogen 
Emission 



Comp aring our results with the Ha distribution of lConselice et alJ 
(2001) (reproduced partly here in Fig.0, we find the molecular 
hydrogen emission is clearly associated with the atomic emission, 
suggesting a common excitation mechanism. 

Molecular to atomic line ratios are commonly used to distin- 
guish between various excitation mechanisms. 3o" upper limits of 
Bry were taken from our spectra and estimates of the Ha emission 
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Region 


Excitation 
Temperature (K) 


ortho-to-para 
ratio \i) 


H2 Column Density 
cm j 


Horseshoe knot 


2050*38 


2.8±fj 


7.7 ±2.0 


SW1 knot 


1720±!1> 


3.0±f;l 


10.1 ±1.6 


Region B 


1660±lfg 


1 9± 08 


13.5±3.1 


Region C 


zlJU=t 440 




2.5 ±1.0 



Table 3. Excitation temperature, ortho-to-para ratios and total line-emitting 
molecular hydrogen column density of all regions observed which revealed 
more than two molecular hydrogen lines. Errors correspond to the lo" level. 



1.0 



0.1 



Region A 



4 6 8 10 12 

v=l-0 S(l) Intensity (10 erg s^cm^arcsec" 2 ) 



Figure 7. The line intensity ratio v=l— S( 1 )/Bry plotted against v=l— 
S(l) intensity. u=l— S(l) line has lo" errors, whilst Br/ line is a 3o" upper 
limit. Region A is labelled. Active galaxies and starforming regions have a 
ratio ^ 1 (below the dashed line). The weakness of the Br/ line rules out 
stellar UV as a possible excitation mechanism 



from Fig. 5 of lConselice et alj 1200 ll) are used in conjunction with 
the strong v=l— S(l) line to constrain excitation models. 

Active galaxies (starburst, Seyfert and ultra luminous IR 
galaxies) and star forming regions have the ratio v= 1 — S( 1 )/Bry ^ 
1 IJaffe et alJ20oTh whereas most of the filament regions observed 
here have v=l— S(l)/Bry ^ 1 (Fig. 0, which strongly argues 
against a stellar UV excitation model. The exception is region 
A, which has a spectrum with v=l— S(l)/Bry^l, although the 
molecular and Bry emission may have different origins as Fig|4| 
suggests. All regions observed, except regi on A, do not show ev- 
idence for stars in the continuum images IConselice etaill200ll) 
and regions further from the galaxy (SW1, SW2 and the horseshoe 
knot) do not show a continuum in their K-band spectra (Fig.[2j. It 
therefore seems unlikely that these regions can be excited by UV 
photons from young, hot stellar populations. However Lya emis- 
sion and UV continuum ha s been found associated wi th the opti- 
cal filaments (Fabian et al. 1984; Norgaard-Nielsen et al. 1990) that 
may not have originated from stars. 

The v=l— S(l)/Ha ratios for the outer filaments, shown in 
Fig- 111 generally lie between 0.02-0.1, similar to ratios from the 
centre of other domin ant cluster galaxies. In the central region of 
NGC 1275 iDonahue et alj i200rl show the v=l-0 S(l)/Ha ratio 
decrea se from 0.19 t o below 0.01 out to a radius of 4 arcsec. Fol- 
lowing IShields et alj ll990h estimate of the ratio of Ha emission 
in the filament compared to the stellar cluster of region A, only 



s 

X 



0.100 



0.001 



bRegion A 



10 20 
Distance from NGC 1275 (kpc) 



Figure 8. The line intensity ratio v=l— S(l)/Ha as a function of distance 
from the nucleus of NGC 1275. The region bounded by the dashed lines 
corresponds to the range o f v=l— S (l)/Ha ratios in other nuclear regions 
of central cluster galaxies IWilman et at]|2002l) . As there is no correlation 
with radius, the central AGN can not be exciting the molecular hydrogen 
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Figure 9. The line intensity ratio Pa/Ha for regions in which Pa could be 
measured or a reasonable estimate made. All errors follow from 1 o" errors 
on the line intensities. The dotted lines give the expected ratio given Case 
A recombination with gas at the temperature written above the line. The 
dashed lines give the expected ratio given Case B recombination with gas 
at the temperature written below the line lOsterbrock 1989). Pa emission 
from region A is divided into the broad and narrow component, with 30 per 
cent of the detected Ha emission assumed to have come from the filamen- 
tary material, A(fil), and 70 per cent from the Stellar Cluster, A(SC). 



30 per cent of the measured Ha surface brightness has been taken 
into account for region A whereas all the v=l —0 S(l) emission has 
been taken into account. Even with this estimate, the ratio of region 
A falls below the other regions. Fig. [8] shows no correlation with 
distance from the central galaxy as would be expected if the active 
nucleus of NGC 1275 enhanced the Hll or H2 emission. We infer 
that the AGN can only pla y a minor role in exciting the filaments 
Jjohnstone & Fabianll98l) . 

Possible excitation mechanisms that could pr oduce these ra- 
tios h ave been discussed in a number of papers. IDonahue et alj 
(2000) suggests X-ray heating of thick columns of molecular hy- 
drogen, conduction from the ICM or mixing layers. These ra- 
tios could also be obtained from shocks according to the models 



© 0000 RAS, MNRAS 000, 000-000 



8 N. A. Hatch, C. S. Crawford, A. C. Fabian & R. M. Johnstone 



of iHollenbach & McKed <1989l) . whilst predictions of CLOUDY 
models, of dense self-gravitating clouds with rij = 10 5 cm~ 3 ex- 
posed to UV radiation also produce these ratios IWilman et alJ 

The Ha intensity of the H? regions was e stimated using aper- 
ture photometry on the lConselice et all 1200 ll) image (Fig0 over 
the regions were the slits crossed. In doing so errors may be intro- 
duced by, for example, picking slightly different regions and differ- 
ent calibration methods. In addition no correction has been made 
to remove the [Nil] line flux from the image or take into consider- 
ation any reddening which could have a large effect on this ratio. 
As it is difficult to assign an error to this, we can instead measure 
the Pa/Ha ratios and compare them with predicted recombination 
ratios. Fig [5] shows the Pa/Ha ratios for regions SW 1, A, B and 
C toge ther with various recombination ratios given bv lOsterbrocld 
1 1989). Although there is a large scatter, regions SW1, B and C all 
agree within error to Case B recombination ratios. The Pa/Ha ra- 
tio for both the broad filamentary emission and the narrow stellar 
cluster emission of region A lies well below any expected recom- 
bination ratio. Therefore it is possible that the v=l— S(l)/Ha ra- 
tio of region A is badly compromised by an inaccurate Ha surface 
brightness estimate as well as an uncertain v=l— S(l) line surface 
brightness. 



4.4 Pressure Balance 

The X-ray emitting gas surrounding the molecular hydrogen has 
a number density of n~0.06cm~ 3 , and a temperature of approxi- 
mately 3x 10 7 K, giving a thermal pressure of ~ 1.8 x 10 6 cm -3 K 
JSanders et all2 004). The pressure of the outer regions of the op- 
tical line-emitting gas has been estimated from the [Sll] doublet 
I Heckman et al. 1989) and the gas is found to be approximately 
in pressure equilibrium with the surrounding X-ray medium. As 
the molecular hydrogen line ratios imply the emission is partially 
if not predominately thermal, with T~ 1 600-2200 K, the total den- 
sity from at least part of the molecular line-emitting region must 
be greater than the critical density above which thermal emission 
dominates, n^l0 4 cm~ 3 . This implies the molecular line-emitting 
regions must be over-pressurized by at least 10 times that of its 
surrounding m edium . This echo es the press ure problem raised by 
Wilm an et alJ <2002h and Ijaffe et al] 1200 ll) who postulated the 
need for a dynamic, non-isobaric model in order to explain both 
the ionized and molecular emission. 



5 CONCLUSIONS 

We have found molecular hydrogen emission from the same re- 
gions as the Ha emission in the outer filaments, up to 24kpc 
from the nucleus of the central cluster galaxy, NGC 1275. As 
H2 quickly dissociates above temperatures of 4000 K, this finding 
proves that a cooler component is intimately associated with the ex- 
tended optical/UV line-emitting filaments. The spectra are consis- 
tent with predominately thermal emission, whilst the ortho-to-para 
ratios of regions B and C suggest an additional non-thermal exci- 
tation component possibly from UV photons in the Lyman-Werner 
band. The high H2/Bry ratios suggest that UV radiation from hot, 
young stars is unlikely, however there is Lya emission and UV 
continuum associated with the optical filaments and X-ray irradi- 
ation of the molecular hydrogen by the surrounding ICM may be 
able generate enough photons in Lya and in the Lyman-Werner 



band through non-thermal electron collisions with hydrogen atoms 
jMalonevetalll99d) . 

The measured excitation temperatures lie between 1600- 
2200 K and the column densities of the warm excited molecular 
hydrogen in these regions are 2— 14x 10 15 cm~ 2 . It is possible that 
there exists an even cooler component of molecular gas within 
these filaments which will be revealed only by searches for the 
v = (0 — 0) pure rotational lines of hydrogen in the mid-IR, or CO 
emission in the mm waveband. As the emission is thermal, at least 
some of the emitting gas has a density greater than 10 4 cm~ 3 . Com- 
paring the minimum pressure of the molecular line-emitting gas 
to that of the X-ray gas shows the molecular gas must be over- 
pressurized by at least an order of magnitude. 

The presence of this molecular gas so far out in the ICM poses 
similar challenges as the presence of the optical line-emitting gas. 
Either the gas has formed in situ (having cooled from the ICM and 
is falling onto the galaxy in a condensing phase), or the gas has 
travelled from a much cooler region, abundant in molecular hydro- 
gen, to where it has been detected in the ICM. 

A large reservoir of molecular hydr ogen is known to lie 
in the nucleus of NGC 1275 ijaffe et al]|200lt lEdee et aljEool 
Donahue et al. 2000;: lKrabbeetall2 000). The origin of this gas is 
still debate d, but it may have c ome from a previous merger as sug- 
gested by Braine et al. 1 1995) or be accumulated through cooling 
of the ICM. Theories have been proposed suggesting that the buoy- 
ant old radio bubbles, visible in X-ray ima ges, drag gas from the 
central regions out into the ICM as they ri se I Bohringe r et aljl995t 
Chura zov et alj |2001 : Fabian et al. 2003). In these models the fil- 
aments act as streamlines tracing the path of the cooler gas as it 
moves and is heated by the ICM. 

The horseshoe-shaped filament and its mirror image seen to 
the West (see Fig. in the NW, appears just below an old radio 
bubble with a flow pattern similar to t hat of water underneath a ris- 
ing air bubble (see lFabian et all2"0 03 for discussion). It is possible 
that the molecular hydrogen originated from the core of NGC 1275 
and it is currently being heated as it is dragged into the hot ICM be- 
hind these buoyant bubbles. However not all optical emission-line 
filaments with molecular hydrogen appear to follow visible radio 
bubbles; in particular the radial filaments and molecular knots in 
the East do not point toward any known radio bubble. These fila- 
ments could have been drawn out by a bubble whose radio emis- 
sion is now too faint to detect, although that would suggest that the 
molecular hydrogen exists in the hostile environment of the ICM 
for about ~ 10 7 — 10 s years, thus having clear implications on the 
rate of conduction within the ICM. High velocity resolution spec- 
troscopy of these filaments would enable mapping of any flow pat- 
tern and thus probe further the origin of this cool gas. 
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